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Inspired by the recent T2K indication of a relatively large #13, we provide a systematic study of 
some general modifications to three mostly discussed neutrino mixing patterns, i.e., tri-bimaximal, 
bimaximal and democratic mixing matrices. The correlation between #13 and two large mixing 
angles are provided according to each modifications. The phenomenological predictions of 612 and 
823 are also discussed. After the exclusion of several minimal modifications, we still have reasonable 
Q,' predictions of three mixing angles in 3a level for other scenarios. 
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00 ' I- INTRODUCTION 
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The observation of neutrino oscillations has revealed that neutrinos may have non-zero masses and lepton flavors 
are mixedflj. In the basis where the flavor eigenstate of three charged leptons are identical to their mass eigenstates, 
the mixing of neutrino flavors can be described by the standard parameterization [2], which is expressed by three 
mixing angles 9 12 , 9 23 , 13 and three CP-violating phases 6,p,a: 



V 



c 12 c 13 s 12 c 13 s 13 e 

— s 12 c 23 — c 12 s 13 S 23 e c 12 c 23 — s 12 s 13 S 23 e C 13 S 23 

S 12 s 23 — C 12 s 13 C 23 e — c 12 s 23 — s 12 s 13 C 23 e C 13 C 23 / 
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where Cy = cosd i ,, s^ = sin^- and P v = {e lp , e ta , 1} is a diagonal phase matrix, which is physically relevant if 
neutrinos are Majorana particles. It is possible to fit neutrino masses and mixing parameters from various neutrino 
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oscillation experiments. From|3j, we have 



a 



Am2! = 7.59±0.20(to;g9) x l(T 5 eV 2 



3^t2 



Am|i = -2.36 ±0.11(±0.37) [+2.46 ± 0.12(±0.37)] x lO^eV 

9 12 = 34.5 ± 1.0(tli)° , ^23 = 42.8^(^° 3 7 ) , 0i3 = 5.1±1;°(< 12.0)° , 



(2) 



at lcr(3cr) level. 

However, the latest result from long baseline neutrino oscillation experiment T2K indicates that 13 is relatively 
large. At the 90% confidence level, the T2K datajj) are 



5-0° < 0i3 < 16.0° ( NH ), 
(IH), 



5-8° < 0i3 < 17 



(3) 



for a vanishing Dirac CP-violating phase 8 = 0°, where "NH" and "IH" correspond to the normal and inverted 
neutrino mass hierarchies respectively. The best fit values are 6,-, = 9.7° (NH) and 11.0° (IH). There are already 
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some literatures discussing on this issue [5J, |6J, I10h16| . 

From the theoretical point of view, there are three types of well motivated neutrino mixing patterns: tri-bimaximal 
mixing pattern (TB) [8], bimaximal mixing pattern (BM)[7j, and democratic mixing pattern (DC) [9j, which may 
arise from some discrete flavor symmetries, such as A 4 and fi — r symmetry, or some very special structures of neutrino 
mass matrices. The explicit forms of them read as follows: 
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It is obvious that all these three mixing patterns are not consistent with the current neutrino mixing data, because 
they all have predictions of vanishing 13 . Actually, the paper [5j considered possible perturbations to the Democratic 
neutrino mixing pattern to get the relatively large 13 . While the paper 6j considered minimal modifications to the 
tri-bimaximal mixing pattern to fit the T2K data. 

In this paper, we consider some general modifications to these neutrino mixing patterns to get appropriate neutrino 
mixing angles that may fit the T2K result. Since C^pmns = UTU vy where U e and U v are used to diagonalize the 
charged lepton mass matrix and neutrino mass matrix respectively, these modifications may come from the neutrino 
sector or from the charged lepton sector, and even both. We will discuss these cases in detail. 

The remaining part of this paper is arranged as follows: In section II, possible modifications are listed. Section III 
is devoted to study their phenomenological results. We summarize in section IV. 



II. POSSIBLE MODIFICATIONS FOR V TB , V BM AND V oc 



In this section, we consider possible modifications for V TB , V BM and T^ DC . For a given neutrino mixing matrix U, 
there are three possible forms of modifications: U ■ X, Y ■ U or Y ■ U ■ X, where X and Y denote generic perturbation 



matrix. Notice that the neutrino mixing matrix, which comes from the mismatch between the diagonalizations of 
the neutrino mass matrix and the charged lepton mass matrix, is given by t/pMNS = ^i *u- U ■ X may come from 
perturbations to the original neutrino mass matrix, which is obtained from certain flavor symmetry, Y ■ U may come 
from perturbations to the original charged lepton mass matrix, while Y ■ U ■ X may come form perturbations to the 
both sector. 

From the mathematic point of view. X and Y can be expressed as X(Y) = V 23 V 13 V 12 in general, where V 23 , V 13 
and V 12 are given by 

COS; 

V 12 = - sin 

V ° 

where x, y, z denote rotation angles and S denote possible CP-violating phases. For simplicity, we may turn off one 
or two V^ (ij = 1,2,3). In this way, we may explore easily the physical meanings of these perturbation, see Ref[5|, |6[ 
for illustration. To be consecutive, we first classify possible simple perturbations and list them below: 
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where a = TB,BM or DC and (ij), (kl) = (12), (13), (23) respectively We left Y ■ U ■ X case for future discussion 
due to its tedious calculation. In the next section, we will investigate the phenomenology of these perturbations, by 
indicating their predication on the most mysterious neutrino mixing angle 9 13 . Then, future neutrino oscillation result 
of 9 13 may also exclude or support these modifications. 

III. PHENOMENOLOGY 

We investigate in this section phenomenologies of perturbation scenarios listed in Eqs. EEl We first investigate the 
most simple cases presented Eqs. [5]and[7J The phenomenologies of scenarios V TB ■ V^ and V^ ■ V TB were already 
studied in Ref. p, such that we will only investigate the other two cases. 

For the V BM ■ V^ and V^ ■ V BM case: 

• V BM • V 12 : We have 9 13 — and 9 23 < 45° in this case. It's thus excluded by T2K. 

• V BM ■ V 13 : We have 9 12 > 45° in this case. It's thus excluded. 

• v bm ' v 23 : Given 9 12 C [31.7°, 37.7°], 9 13 lies in the range [26.7°, 33.7°]. It's thus excluded. 



• ^23 ' V BM : 9 13 = and 9 12 = 45°. It's thus excluded. 



• V 13 • V BM : In this case, we have 



tan 8 12 = 
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In the left panel of Fig. IIILll we plot 9 12 , 8 13 and 6 23 as functions of z. For simplicity, we choose 5 = 0. 
that there are enough parameter space for this scenario to fit both the T2K's result on 9 13 and other 
oscillation data. For the case S ^ 0, we can derive the Jarlskog of this scenario. 

sin 2z (2 -sin 2 z) (3 + cos 2z) 2 - 8 sin 2 2z cos 2 S . 
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4^2(1 + cos 2 z) (3 + cos 2z) 2 + 8 sin 2 2z cos 2 8 
Given 9 13 = 9.7°, we plot in the right panel of Fig. IIII.ll the Jarlskog as function of S. We may read 
figure that J lies in the range [0, 0.08] in this scenario. 



(10) 

(11) 
(12) 

It's clear 
neutrino 

(13) 
from the 













V13.VBM 


50 












40 






e 23 






30 












20 




e 12 








10 




i i 


■ 


9l3 

■ 






FIG. III.l: Case BM1. 0%j as a function of 8 (left panel) and Jarlskog invariant as a function of S (right panel). 



V 12 • V BM : We have the following correlations: 



tan 8 12 = 
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cos z + sin ze 
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To generate proper 8 12l we set 6 = it. We may find that this scenario is quite similar to the last one. The only 
difference is that this scenario predict 8 23 < 45°, while the last scenario predict 8 23 > 45°. Both scenarios are 
workable. 

For the V DC ■ VL and V^ ■ V DC case, 



• V DC ■ V 12 : 9 13 — in this case. It's thus excluded. 



• V DC • V 13 : We have # 12 > 45° in this case. It can be excluded. 

• V DC • V 23 : Given 6 12 C [31°, 38°], 9 13 lies in the range [26°, 35°] in this scenario. It's thus excluded. 



• ^23 ' V DC : 13 — 0. It's excluded. 



• V 13 • V DC : # 23 > 54.7° in this scenario. It's thus excluded. 



• V 12 • V DC : In this scenario, we have sin# 13 = y/{2 — tan 2 23 )/3. Given # 23 c [35.5°, 53.5°], 8 13 lies in the 
range [13.5°, 45°]. It's excluded by the best fit value but allowed by T2K's data. If 8 23 is precisely given in the 
future neutrino oscillation experiments, We may get to the denial of this scenario. It would be worth mentioning 
that proper 9 12 can always be obtained in this case. 

Now we go to investigate phenomenologies of a little more complicated cases: V a ■ V^ ■ V kl and V^ ■ V kl ■ V a 
(a = BM,TB,DC). 

• VTBrrl = VtbV 23 Vi 3 : We get the following correlations 
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Here sin y can be expressed as functions of 9 13 and 6 12 
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1 - 2 tan 2 6<i2 + 3 sin 2 6 13 



y 1 + tan^ 9 12 

Then we can get sinz from the relation of #13 with y and z, i.e. 



(20) 



smz = 



V6 sin #13 + sin y J 2 — 3 sin 9\ 3 + sin 2 y 



(21) 



2 + sin 2 y 

Substituting Eq |20l and [5T] into tan #23 , we can obtain the relations between #23 and #13 which is shown in 
Fig |III.2l We can see from Fig. IIII.21 that as #13 increases from 0° to 18°, decreasing #23 from 50° to 39° could 
be get. Both #13 and #23 are predicted in the experimentally allowed region. 

• VTBrr2 = V T bV 2 3Vi2 : We have 
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FIG. III.2: Case VTBrrl. 6*23 as a function of 0i 3 . 




FIG. III. 3: Case VTBrr2. 813 as a function of 823 (left panel) and 812 as a function of x (right panel). 



From the above equations, we could get 



sin 0i3 



1 — tan 1 
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^5(1 + tan 6j 3 ) -2tan0 23 



(23) 



The illustration of #13 as a function of 023 and $12 as a function of Euler angle x is shown in Fig. 1111.31 We can 
see from the left panel that this kind of modification from tribimaximal mixing pattern provides the prediction 
of 0i3 value up to 6° while the 023 varies in the region allowed by current experimental data. In the right panel, 
reasonable 0i2 could also be produced in this case when rotation angle x is changed. 



• VTBrr3 = VTBV13V12 : For this scenario, we have 
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Thus #13 and 6*23 have the following correlation: 
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The numerical results of this scenario is shown in Fig. IIII.4I 
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FIG. III. 4: Case VTBrr3. #13 as a function of #23 (left panel) and #12 as a function of x (right panel). 

In this case, #13 varies from 0° to almost 12° when #23 is set in its experimentally allowed region. In addition, 
when we change the corresponding rotaion angle x from 0° to 10°, #12 varies from 35.7° to 38.2°, which is 
consistent with its global fit data in 3tr level. 



VTB111 = V23V13VTB : In this case we have 
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Then 6*13 could be represented by 6*12, i.e. 
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FIG. III. 5: Case VTB111. #13 as a function of #12 (left panel) and 623 as a function of y (right panel). 



From Fig. IIII.5I we can see in the left panel that #13 is less than 6°in this modification scenario from the 
tribimaximal mixing pattern. In the right panel, #23 from 34° to 46° could be produced which lies in the region 
of global fit data [3| in 3ct level . 



• VTB112 = V 23 V 12 V T b : 
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Thus, 013 could be produced, i.e., 
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'4 + 5 tan 0^ 2 - 2v / 2tan0i 2 

Numerical results are illustrated in Fig. IIII.61 it is obvious that this case shares the same 0i3 prediction with the 
last case, which could only provides very little 0i3, i.e. less than 6°. But 023 evolves differently, the prediction 
of this mixing angle could be 44° at most, which is 2° smaller than the last scenario. 



• VTB113 = Vi 3 V 12 V T b 
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FIG. III. 6: Case VTB112. #13 as a function of #12 (left panel) and #23 as a function of y (right panel). 



cos a: could be deduced from the above equations, i.e., 
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so is the cosz: 
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Numerical results is shown in Fig IIII.71 In this case, as #13 increases from 0° to 18° we can get #12 from 35° 
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FIG. III.7: Case VTB113. 6 12 as a function of 6>i 



to 5°. It is obvious that #13 less than 5° is favored if we choose only reasonable #12 from 30° to 35° in the 3<r 
region. 



For the scenarios V BM ■ V^ ■ V kl and V^ ■ V kl ■ V BM , we have 
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• VBMrrl = V BM ■ V 13 ■ V 12 : For this case, we obtain 
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Since both 9 13 and (9 23 depend on a single parameter y, they have the following correlation: 
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In the left panel of Fig. IIII.81 we plot 9 13 as a function of 9 23 . Taking into account the 3<r constraint on 9 23 , 
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FIG. III. 8: Case BMrrl. 613 as a function of 623 (left panel) and #12 as a function of y (right panel). 

the 6 13 predicted by this kind of modification lies in the range [0°, 9.5°], the upper bound of which is approach 
to the T2K's best fit value for normal hierarchy case. Appropriate value of 9 l2 may be obtained by varying y, 
which is shown in the right panel of Fig. IIII.8I 



• VBMrr2 = V BM • V 23 • V 12 : In this case, we have 
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It's easy to check that the correlation between 6* 13 and 6 23 is the same as that of Eq. [23 The only difference 
between this scenario and the last one is that their 9 l2 evolve differently. Both scenarios predict reasonable 9 12 
with a big rotation angle x. 



• VBMrr3 = V BM • V 23 ■ V 13 : For this case, we derive the following a little complicated correlations 
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■ is a little more complicated. Here, we only show their numerical results. 



11 

(42) 
(43) 

(44) 

(45) 

(46) 

In the 



50 



45 
40 
35 
30 
9 " 25 
20 
15 
10 



VbM-V23-Vi3 



_l . I . I . I . I . I . I . I . I . 1_ 



10 12 14 16 18 20 



FIG. III.9: Case VBMrr3. 6» 23 as a function of 6» 13 



left panel of Fig. IIII.91 we plot 9 2 s as function of 9 13 , choosing the best fit value of 9 12 . We may read from the 
figure that the upper bound for the 9 13 is 6.1° for this case. 



• VBMftl = V 13 • V 12 • V BM : It has the following correlations: 
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Here x and z can be expressed as function of Q i 
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We show in Fig. IIII.101 9 12 as function of 8 13 , setting 9 23 its best fit value. We may find that 9 13 lies in the 



45 




Fi3.Fl2.Vbm 


40 






35 






30 






25 






20 






15 


. i . i . i . i 


. i . i . i . i . i . i . 



10 12 14 16 18 



20 



FIG. III. 10: Case BM111. 8 12 as a function of 0i 3 . 
range [3.7°, 7.2°], given 9 12 C [31.5°, 37.5°]. It's impossible to generate larger 9 13 in this scenario. 
• VBMI<!2 = V 23 • V 13 • V BM : We may derive the following relationships from this scenario 
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It's easy to check that 8 13 and 9 12 have the following correlation: 
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Since the value 9 l3 only depend on 9 12 , we may constraint the 9 13 s parameter space, according to this equation. 



It's numerical result is shown in Fig. MI. Ill We may read from this figure that 9 13 lies in the range [6.4°, 10.6°]. 
We plot in the right panel of Fig. 1111.111 9 23 as function of y. Appropriate 9 23 can be obtained by small 
perturbation to y. 
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FIG. III. 11: Case BM112. 9n as a function of #12 (left panel) and #23 as a function of y (right panel). 



• VBM^3 = V 23 • V 12 • V BM : The following equations arise in this scenario 
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where 8 12 and 8 13 has the following correlation: 
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As shown in Fig |III.12l #13 decreases with the increase of #12 in this case. When 812 is set to lie in its 3er global 




FIG. III. 12: Case VBM113. #13 as a function of 812 (left panel) and 823 as a function of Euler angle y (right panel). 

fit region, 813 ranging from 8° to 15° could be predicted which lies in the larger part of the T2K allowed region. 
When Euler angle y varies, 823 changes from 44° to 54° which is consistent with its 3a global fit data. 



For the scenarios V BM ■ V i j ■ V kl and V^ ■ V kl ■ V BM , we have 
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• VDCrrl = V DC • V 13 • V 12 : In this scenario, we have 
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where 9 13 and # 23 have the following correlation: 
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We plot in the left panel of Fig. 1111.131 9 13 as function of 9 23 . In this scenario the possible parameter range 





FIG. III. 13: Case VDCrrl. 9i3 as a function of 822 (left panel) and #12 as a function of y (right panel). 

of 9 13 is [3.5°, 16.5°]. Big 9 13 is permitted in this scenario. If future reactor neutrino oscillation experiment 
constraint 9 13 lying below 3°, this scenario can be excluded. We also plot 9 12 as function of x, given 9 it's best 
fit value 9.7°. We may obtain appropriate 9 with big rotation angle x. 

• VDCrr2 = V DC • V 23 • V 12 : The phenomenology of this scenario was already studied in paper [5j. We won't 
repeat them here. 



• VDCrr3 = V DC • V 23 • V 13 : We may derive the following equations in this scenario 
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Given these correlations, we may plot 9 23 as function of 9 13 by choosing 9 12 its best fit value. From Fig. IIII.141 
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FIG. 111.14: Case VDCrr3. 9 23 as a function of 13 . 

we may read parameter space of 9 13 by setting 9 23 changing in it's 3c range, which is [0.9°, 11.1°]. It's an 
excellent parameter space, covering T2K's best fit value for both normal hierarchy case and inverted hierarchy 
case. 

VDCffl = V 13 • V 12 • V DC : In this scenario, we have 

V6 cos x cos z + V2 sin x cos z + 2 sin z 



tan# 12 



sin I 



13 



y/6 cos x cos z — v2 sin x cos z — 2 sin z 
%/2 sin x cos z — sin z 



tan 9. 



23 



V^cosx 



cos z + V2 sin x sin z 
Here, x and y may be expressed as functions of 9 {j 



SIM 



sinz 



'2cot 2 2 3 + 3sin 2 13 - 1 

2(1 + cot 2 9 23 ) ' 

tan 2 3 (tan 12 — 1) — sin0 13 (l H-tan^ 12 ) 



(68) 
(69) 

(70) 



(71) 
(72) 



3(1 + tan6* 12 ) + yosin^g tan0 23 (l — tan(9 12 ) 
We plot in Fig. IIII.15I 2 3 as function of 9 13 , by setting 9 23 its best fit value. It's clear that we can not get the 
appropriate 9 12 when 9 13 changes in the range [0°, 18°] in this case. Such that it's excluded. However, if we 
let 9 23 changing in its 3er range, a narrow appropriate parameter space of 9 12 can be obtained. In short, this 
scenario doesn't work for the best value. 
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FIG. III. 15: Case VDC111. 13 as a function of 8 12 . 



VDCII2 = V 23 • V 12 • V DC : In this case we have 



tan0 15 



smf 13 = 



a/3 cos X 


+ sin 2 


a/3 cos X 

\ll s[nx 


— sinx 



tan(9. 



23 



a/2 cos a; cos y + sin y 



cosy 



V2 



cos a; sm y 



The correlation between 13 and 12 is 

sin 6*i 



1 — tan ( 



(73) 
(74) 
(75) 

(76) 



a/2 + 2 tan 2 12 - 2 tan 6» 12 

We plot in the left panel of Fig. IIII.16[ 13 as function of 12 and in the right panel of Fig. IIII.16I 2 3 as function 
of y. We find that 13 can only changes only in the range [12.5°, 21.3°]. This range is excluded by the best 
fit value, but allowed by the T2K's result. If future reactor neutrino oscillation experiments negate the T2K's 
result, this scenario will be excluded definitely. 



• VDCM3 — V 23 • V 13 • V DC : The following equations can be obtained from this scenario 

a/3 cos z + a/2 sin z 



tan0 12 



sm( 



13 



tan0, 



a/3 cos z — a/2 sin z 
smz 

a/2 cos y + cos z sin y 



cos y cos z 



V2 



sin y 



We have 



1 — tan0 l: 



a/5 tan 2 12 — 2 tan 12 + 5 



(77) 
(78) 
(79) 

(80) 
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We plot in Fig |III.17l fl 13 as function of 9 12 and 9 23 as function of y. It's clear from the figure that 9 13 can only 
change in the space [6.4°, 10.5°] in this scenario. We can get appropriate 9 23 with big rotation angle y, which 
can be the result of certain flavor symmetry. 
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FIG. III. 16: Case VDC112. #13 as a function of #12 (left panel) and #23 as a function of y (right panel). 




FIG. III. 17: Case VDC113. #13 as a function of #12 (left panel) and O23 as a function of y (right panel). 



IV. SUMMARY 



Among the knowns and unknowns of neutrino physics, the nonzero smallest lcpton mixing angle #13 has received 
a lot of attentions. The studying and measurement of #13 is important to enrich our understanding of neutrino 
properties. Inspired by the recent T2K result of a relatively large 9 13 , We investigated some general modifications, 
which may minimal modifications or non-minimal modifications to the three well-known neutrino mixing pattern: 
Bimaximal, Tri-bimaximal and Democratic mixing patterns. Some non-trivial correlations between 9 12 , 9 23 and 9 13 
were obtained. We investigated constraints on the 9 13 by these correlations. Some modifications are already excluded 
by the current neutrino oscillation data, while the others give their predications on 9 13 . Future neutrino oscillation 
experiments may confirm or negate these models. Since 9 13 is not so small, investigating leptonic CP-violating effects 
become relevant, but beyond the scope of this preliminary work. We will study this interesting and important topic 
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somewhere else. 
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